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Abstract 

The 7-particle form factors of the fundamental spin field of the 0(3) 
nonlinear a-model are constructed. We calculate the corresponding 
contribution to the spin-spin correlation function, and compare with 
predictions from the spectral density scaling hypothesis. The resulting 
approximation to the spin-spin correlation function agrees well with that 
computed in renormalized (asymptotically free) perturbation theory in 
the expected energy range. Further we observe simple lower and upper 
bounds for the sum of the absolute square of the form factors which may 
be of use for analytic estimates. 



1 Introduction 



1 + 1 dimensional integrable models are theoretical laboratories where one can study 
selected problems in field theory. These include the question whether the continuum 
limit of the lattice regularized version of a field theory coincides with the theory 
constructed by continuum techniques. One can study here structural questions 
related to asymptotic freedom, topological excitations etc., and, in the case that the 
theories agree, the nature of lattice artifacts. 

The form factor bootstrap method Ej of constructive field theory has been 
used to construct many integrable models. The method consists of three main 
steps. The starting point is the exact scattering matrix that is calculated (sometimes 
conjectured) by the bootstrap method jH]. Next the form factors (matrix elements 
of local operators between scattering states) are calculated by solving the form 
factor (FF) axioms 0. Finally correlation functions are constructed by inserting a 
complete set of intermediate states between the local operators. 

The program has so far only been carried out completely for the Ising model 
j3]. The S-matrix is known for many integrable models and if there are no internal 
degrees of freedom also the many-particle form factors are easily constructed as a 
multiple product of the basic two-particle form factor (which is known again for 
most models). On the other hand, for models with internal degrees of freedom the 
solution of the FF axioms, matrix functional difference equations in this case, is 
not known in general. Nevertheless, the form factors of some models belonging to 
this class, the Sine-Gordon model, the chiral Gross-Neveu model, and the 0(3) and 
0(4) nonlinear cr-models, are available [21 E]- Finally (with the exception of the 
Ising model) the r-particle contributions to the correlation functions can only be 
computed numerically and therefore this sum has to be truncated after the first few 
terms. 

Among models with internal (isospin) degrees of freedom, it is the 0(3) model, 
where one can go furthest. This is because in this special case the many-particle form 
factors (after removing a simple factor) are polynomials in the particle rapidities 
which can be determined recursively [HI. Fortunately it is also a model which is 
of particular interest since it has many properties akin to Yang-Mills theory in 4 
dimensions in that it exhibits asymptotic freedom and has instanton solutions. 

The question mentioned above concerning the equivalence between the FF con- 
struction and that from the lattice regularization has been addressed in ref. jZj . The 
two constructions are both non-perturbative and in both cases it is the low energy 
properties that are most easily accessible. However in the bootstrap approach it is 
simpler to obtain reliable results at higher energies. For example for the Fourier 
transform of the 2-point correlation function of Noether currents, the contributions 
of higher r-particle states become significant at typical energies rapidly increas- 
ing with r. Of course, despite the fact that in practice renormalized perturbation 
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theoretic computations often seem to be quite accurate down to unexpectedly low 
energies, the property of asymptotic freedom is a statement pertaining to asymp- 
totically high energies, and to establish this property in the bootstrap framework 
it is necessary to have rigorous knowledge on the contributions of r-particle states 
with arbitrarily large r. Although the latter has not yet been achieved, based on 
the exact form factors up to 6 particles. Max Niedermaier and one of the present 
authors (J.B.) presented convincing evidence for remarkable scaling properties of 
the higher intermediate particle contributions jHl IH], which are very powerful be- 
cause they make it possible to include dominant contributions from all states in the 
form factor expansion of correlation functions. Furthermore the scaling hypothesis 
is completely compatible with asymptotic freedom. 

Although the lattice regularization is an important non-perturbative method 
(and practically the only one available in 4d QCD), there are many physical phe- 
nomena which are inaccessible in this approach for example nuclear structure func- 
tions at small Bjorken x. In the framework of the bootstrap approach in 2d however, 
we hope that properties of the structure functions at small x can be extracted, an 
insight which would hopefully also be relevant for QCD. This hope would be realized 
if we could find evidence for scaling properties of the structure functions similar to 
those already observed for the 2-point function spectral densities. 

With this goal we have initiated a project to compute the structure functions 
in the 0(3) sigma model, and with the known form factors we have computed 
contributions from intermediate states up to 5 particles. To establish their scaling 
properties however it is helpful to extend the list of known form factors to higher 
particle number. In this paper we discuss the calculation of the 7-particle form 
factors, which although in principle trivial, is technically challenging because for 
r = 7 one has to deal with quite large polynomials. We will report on the structure 
functions in a future publication. Here we restrict attention to the contribution of 
the 7-particle intermediate states to the spin-spin correlation function, which yields 
an additional test of the scaling hypothesis. We also exhibit simple bounds on the 
square of the known form factor polynomials, which if they could be generalized to 
an arbitrary number of particles might be useful for analytic estimates to establish 
general properties. 



2 Some basic definitions 

The central object of attention in this paper is the 7-particle form factor of the spin 
operator of the 0(3) model. This is the r = 7 case of 

(0|5'^(0)|ai, 9,;...; a„ 9,.) = ^ (^i, ...,9^), (1) 
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where the particles are labelled by their isospin indices and rapidities and the spin 
operator S*" also carries an isospin index. The normalization factor 2/1/7? is intro- 
duced here for later convenience. The spin operator itself is normalized by 

(0|5'^(0) |ai, e,) = , (6, x|a, 6) = An6^'6ix - 9) , (2) 

which means that the Minkowski propagator has unit residue at = M^, where M 
is the physical mass of the 0(3) particle. 
Introducing the squared form factors 

5'^'F^^\U) = J2 fl-aAdu---,Orrfl..a,.{0U---,Or) (3) 

the spectral density is given by 

. oo 

p^P^°(^) = -^p(2^+i)(/i), (4) 

TT ^ — ^ 

k=0 

where 

p«(p) = r 5(p - M«(.)) . (5) 

Jo (47r)' 

Here Uj = 9j — 9j^i are rapidity differences and M^^\u) is the r-particle invariant 
mass. The Fourier transform of the correlation function is represented as the Stieltjes 
transform of the spectral density: 



/--(p^) = / dp . (6) 

+ p^ 



We will compare (the truncation of) (jHl) to the results of perturbation theory. 
The 2-loop order perturbative result is |10j 



p2rP-(p2) = Ai|-^ + (2 + ^o) + (2 + eo)«(p) + 0(a2(p)) |> . (7) 
Here the running coupling function a{p) is the solution of 

1 1 / N 1 P 

+ m a [p) = m ■ 



a{p) ' ' M ' ' 

and the parameter gives the connection between the perturbative mass parameter 
and the exact mass gap M. In the 0(3) model their ratio is known exactly [TT] : 

M 

Co = In -— = In 8 - 1 ^ 1.07944 . (9) 
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The overall constant Ai cannot be calculated in perturbation theory, but has been 
determined exactly |H] using the scaling hypothesis (see also Sect. 5): 



A. = (10) 

Instead of the physical form factor (P) it is convenient to consider the reduced 
form factors Qa^ a^ defined by 

...A) = ^^(^1, . . . , e.) gl...,M, ...,er). (11) 

Here 

■^{9^,...,er) = llm-o,), (12) 

i<j 

with 

^W = ,,^"" tanh^^ (13) 

We note that (fTTIjl is the basic 2-particle form factor and the product (fT^ would be 
the complete r-particle form factor if there were no isospin degrees of freedom. 

The advantage of using the reduced form factors is that they are polynomials 
and can (in principle) be recursively calculated The simplest nontrivial case is 
r = 3: 

Ca.aM = <5'^'^M'^l'^2 (02 - 0,) + ^'^-^-1-3 _ Q,^ _ 3^7^) + J^lJ^^aa _ Q^y (14) 

The higher reduced form factors very rapidly become complicated because of the 
large number of isospin components and rapidity variables. 

The quantity entering the spectral densities and two-point functions is the ab- 
solute square of the form factors, summed over the internal symmetry indices. For 
the reduced form factors the corresponding quantities are 

G^^\e,,...,er) = \ K,-aSd,,---A)?. (is) 

For the r = 3 example we get 

^(=^^^1,^2,^3) = 2 [(^1 - e^f + {6, - e^f + (02 - e^f] + i2t:\ (i6) 
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3 The 7— particle form factors 



The 7-particle reduced form factors 



9, 



ai ...ay 



(^l,...,^7) 



(17) 



are polynomials of degree 15 in the rapidity differences. To count the independent 
isospin components we can obviously fix the operator index a = 3 since the rest of 
the components are related by simple isospin transformations. Then we get 274 non- 
vanishing components not counting terms related by the 1^2 isospin symmetry 
twice. These components can be calculated from the known 5-particle form factors 
using the inhomogeneous FF axioms |6j . The resulting polynomials contain typically 
~ 10^ terms. Actually it is not necessary to do this calculation for all the 274 
components since we can find 4 components such that the rest can be obtained 
from them by using the homogeneous FF axioms (corresponding to permutations 
of the rapidity variables). After having computed all terms we checked once more 
that all the FF axioms are satisfied by p7|) and also that the coefficients of the 
leading terms in the last variable (those proportional to 9^) are proportional to the 
6-particle reduced form factors as they should jH] . 

After the calculation of the form factors the next step is to compute the sum 
of the absolute square of the components since this is needed in the calculation of 
the 7-particle contribution to the spin-spin correlation function and other physical 
quantities. This calculation is rather challenging because at intermediate steps one 
has to handle quite large polynomials (of degree 30 in 7 variables consisting of 
about ~ 4 ■ 10^ terms). The final result can be simplified enormously because, 
being a polynomial symmetric under permutations of the rapidity variables, it can 
be expanded in terms of the basic symmetric polynomials cxi, (T2, . . . , o"7 (defined in 
Eq. (A. 3) of ref. |HI). Written in this way, the final result consists of only 3214 terms 
and is reduced to a manageable size. ^ 

To be sure that we have obtained the correct result for the square G^'^\9i, . . . ,67) 
we have performed a number of checks. First we checked that the leading terms in 
the last variable (the coefficient of ^7°) are twice the analogous 6-particle square 
G^^\9i, . . . ,9q) |H]. Next we checked that the overall leading terms (the terms of 
total degree 30) agree with those of a simple symmetric polynomial P^'^\t) [8J. This 
is defined by (for r particles) 




(18) 



perms 



where 



pi>'\t)= n m-9,r+tn^]. 



(19) 



j~i>i 



Which is however stiU too large to usefully reproduce here. 
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In (fTHjl the sum extends over the r! permutations of the rapidity variables 6i, . . . ,6r 
and in ()19|) the correction terms depending on the parameter t do not affect the 
overall leading terms, they are included here for later convenience. Finally, we 
verified that G^'^\9i, . . . ,67) vanishes at the points [Hj 

- 02 = 02 - = iir, Oj = arbitrary j > 3. (20) 

The fact that G'-^-* passed all these nontrivial tests gives us some confidence in 
the correctness of our results. 

4 The spin— spin correlation function 

Having calculated the square of the 7-particle form factors the corresponding con- 
tributions to the spin-spin correlation function and the spectral density can now 
be computed straightforwardly. The only difficulty is that to get correct numerical 
results one has to use high precision arithmetic in order to avoid huge rounding 
errors in calculating G^'^\ As explained above, to put it into a manageable form, 
we expressed it in terms of the basic symmetric polynomials ai, . . . , cr7. Although 
it is a sum of absolute squares and thus it is obviously positive, it is not manifestly 
positive in this form. Actually we experienced large cancellations between positive 
and negative contributions. After rescaling all rapidities by vr all coefficients of the 
polynomial become integers and we can illustrate the large degree of cancellation by 
considering the following integer rapidities: 

61 = 12, 62 = 11, ^3 = 10, ^4 = 9, e, = 8, ^6 = 7, 07 = 0. (21) 

We can exactly calculate the sum of positive and negative contributions in this case 
and we get 

positive part : 26512984286926120356867283491794696305689824 
negative part : 26512984286926120356855853186993819121689824 
total : 11430304800877184000000 

showing a 21-digit cancellation here! 

We have computed the 7-particle contribution to the spin-spin correlation func- 
tion using the VEGAS integration routine, and a subroutine computing G*^''^ invoking 
quartic precision (32 digit) arithmetic. The results, compared with the prediction 
of two-loop perturbation theory, are shown in Figure H Note that here we have no 
free parameters at our disposal in the perturbative calculation as would be the case 
in most other models. As discussed above, we know the exact relation between the 
perturbative A parameter and the particle mass M and also the absolute normal- 
ization of the perturbative curve is available in this model. The form factor results 
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Figure 1: Plot of p^/''P™(p^) against ln(p/Af), showing a comparison of the form factor approach 
to 2-loop perturbation theory. The solid curve is the FF approximation to p^PP^'^{p^) obtained 
by including all intermediate states with < 7 particles. The lower curves show the analogous 
approximations with < 5, < 3,1 particles respectively. The normalization of the (dashed) PT 
curve is fixed according to ifTIUl . Finally the top curve is the PT curve multiplied by a factor 1.05. 

are in very good agreement with perturbation theory in the expected (high) energy 
range. The small deviation for energies p/M ~ 10^ can be accounted for by the 
contribution of r > 7 intermediate particles. To illustrate the fact that this good 
agreement is quite nontrivial, in Figure Q we also show the perturbative result with 
the overall factor changed arbitrarily by 5%. 

5 Scaling 

The 7-particle results also corroborate the scaling hypothesis |H] for the spectral 
densities. To study this aspect we introduce the modified r-particle spectral density 
depending on the logarithmic variable x by the definition 





(22) 
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These are defined for all r > 2, where the cases with r even refer to the spectral 
densities of the two-point function of the isospin Noether current. The graph of this 
function is a bell-shaped curve starting as zero at a; = Inr, reaching its maximum 
at a; = ^^^^ and then slowly decreasing for larger x. Let us introduce the 
rescaled spectral density Y^"^^ by 




Figure 2: Illustration of the self-similarity property of the rescaled spectral densities. The plots 
show Y^^\z) (dashed) compared with Y^^~^^\z) (solid) for r = 3,4,5,6. 

It has been found (based on the study of up to 6 particles) that the shape of the 
rescaled spectral density and the parameters M^'^^ satisfy self-similarity 

lim Y'-''\z) Y{z) (24) 

with universal shape function Y{z) and asymptotic scaling, 

M^'') ~ M, r-^ , {(''^ ~ r^+° (25) 

for large r, with some coefficients M*, C,* and exponents 7, a. Whereas the properties 
of the form factors are consistent with 7 = 1 only, the other exponent can only be 
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determined numerically with result a = 0.27. For those readers unfamiliar with this 
model we point out the amusing fact that the self-similarity holds for all r > 2 and 
thus interrelates the spectral functions of two different isovector operators, the spin 
field and the conserved vector current. Similarly (25) is equally valid for even and 
odd r values, provided we use the normalization introduced in (1). 

Self-similarity continues to be satisfied very well as demonstrated in Figure |21 
To test asymptotic scaling, we used the fitted numerical values based on results for 
up to 6 particles to "predict" for the r = 7 case 

MiP= 0.03188 and 4^^ = 17.73, (26) 

and compared it to the values directly determined from our 7-particle results 

M(^) = 0.03189 and ^^^^ = 17.77. (27) 

Similarly for the integrals 

poo i-oo poo 1 

c(^)=/ dfip^'\fi)= dxi?(^')(x) /i('^)=/ — i?W(x) (28) 

Jm Jo Jo ^ 

we predict 

c^JJ = 1.464 and /i^? = 0.04879, (29) 

and actually get 

c^^) = 1.46(1) and h^^^ = 0.0488(1). (30) 



6 Ultra— positivity 

Because of the large cancellation between positive and negative terms in the rep- 
resentation of G^^^ in terms of symmetric polynomials, it is natural to ask if there 
exists an alternative representation that is manifestly positive. It is indeed possible 
to find such representations. If we arrange the rapidities in decreasing order (which 
is always possible since the polynomial is symmetric under permutations) then all 
Uj = 9j — 9j+i rapidity differences are positive and we found (for all available particle 
numbers 3 < r < 7) that if we expand G^^'' in terms of these differences then all co- 
efficients are positive. Moreover, it is possible to find upper and lower bounds both 
of which are of the simple form (fTHj) such that, expanded in terms of the rapidity 
differences Uj, each and every term in the expansion of G^''^ has smaller coefficient 
than the corresponding one of the upper bound P'^'^^tu'^) and similarly larger than 
the coefficient of the corresponding term of the lower bound P^'^^t^'^^). 
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r 


1 


2 


3 


4 


5 


6 


7 


C(r) 


0.785 


1.009 


1.140 


1.242 


1.327 


1.400 


1.46 


cW(2) 






1.140 


1.206 


1.229 


1.225 


1.200 


r 


8 


9 


10 


11 








cW(2) 


1.164 


1.118 


1.067 


1.013 









Table 1: Numerical values of the integrals c^''-', c^''-'(2). Numerical errors are estimated as ±1 on 
the last digit quoted. 

The value of the parameter characterizing the lower bound is uniformly t[^^ = 2 

(r) 

for all cases 3 < r < 7, whereas the tu giving the lowest upper bound are 

t(f) = 2, ti'') = 9/4, tf = 2.434, = 2.576, 4^^ = 2.688. (31) 

The last three values are numerically approximate ^ and come from the requirement 
that the constant term of P^'^\tu^) be larger than that of C^*"^. This is clearly 
necessary, but at least in the available cases for r > 5 it is also sufficient ^. 

The integrals over the spin and current spectral densities are given as sums of 
the integrals defined in 



d/xp^P'°(/i) = -^c(2^-+i), 

^ k=0 

/•OO ^ 

C="" = / d/ip^"'^'^(/i) = ^c(2^+2)_ (32) 

k=0 

An outstanding question is whether these are finite or infinite in the FF construction. 
Certainly renormalized perturbation theory predicts that (7*^^™ = oo = C'^'^". Also 
the validity of the scaling hypothesis requires this because in this scenario c^'^^ grows 
as ~ r""; the numerical evidence thereof is shown in Figure 01 

If the bounds we found above for 3 < r < 7 continue to be true for all r, then 
the simple structure of (fTH|) allows to study the structure of the correlation function 
analytically. The existence of an upper limit of simple form may help proving the 
existence of the correlation function whereas the existence of the lower bound may 
facilitate the construction of a proof that C'^p™ = oo = C™"^"^ independent of the 
validity of the scaling hypothesis. 



^e.g. the exact value for the case r = 5 is tu^ = 208^/^ 

^The case r = 4 is an exception since the requirement in this case yields a value (34/3)^/"^ 
2.246, which is not quite sufficient. 
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We are not going to discuss these questions further in this paper, but to illustrate 
the usefulness of the existence of the simple lower bound we have calculated the 
integrals c^^\2). These are defined analogously to the ones in p8|l using 
instead of . The simplicity of the integrand allows us to calculate these integrals 
numerically quite effortlessly ^ up to r = 11. The results of the numerical integration 
are given in Table 1. 

Although the integrals c^'^\2) are (apparently) decreasing with r, it is possible 
that they are decreasing slow enough to make the series (j32|) diverge. Indeed, as 
seen in Figure El rc^^\2) is monotonically increasing for all r evaluated up to now. 
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